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Popular Summary 

Satellite remote sensing is such a complex task that until now it could be done 
only assuming that clouds are homogeneous slabs. This assumption does not consider 
that horizontal cloud variability may influence the amount of sunlight clouds reflect 
toward a satellite. This study examines whether horizontal cloud variability causes 
satellite estimates of cloud optical thickness to yield different results if clouds are viewed 
from different directions. The statistical analysis of water clouds (excluding ice-clouds 
for simplicity) in a one year long global dataset of observations by the Moderate 
Resolution Imaging Spectroradiometer (MODIS) reveals that while the estimated optical 
thickness values are consistent for all view directions if clouds are homogeneous, the 
values are much higher for oblique views than for overhead views if clouds are 
heterogeneous and the sun is fairly oblique. After considering a variety of possible 
scenarios, the paper concludes that the most likely reason for estimating larger optical 
thicknesses at oblique views is the enhanced viewing of cloud sides from oblique 
directions. The results will help understand the uncertainties cloud variability introduces 
into satellite estimations of optical thickness. They complement the uncertainty estimates 
that will start accompanying MODIS cloud products in the near future and may 
eventually help correct for the observed view angle dependent biases. 
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Abstract 


This study examines whether cloud inhomogeneity influences the view angle 
dependence of MODIS cloud optical thickness (x) retrieval results. The degree of cloud 
inhomogeneity is characterized through the local gradient in 11 pm brightness 
temperature. The analysis of liquid phase clouds in a one year long global dataset of 
Collection 4 MODIS data reveals that while optical thickness retrievals give remarkably 
consistent results for all view directions if clouds are homogeneous, they give much 
higher x-values for oblique views than for overhead views if clouds are inhomogeneous 
and the sun is fairly oblique. For solar zenith angles larger than 55°, the mean optical 
thickness retrieved for the most inhomogeneous third of cloudy pixels is more than 30% 
higher for oblique views than for overhead views. After considering a variety of possible 
scenarios, the paper concludes that the most likely reason for the increase lies in three- 
dimensional radiative interactions that are not considered in current, one-dimensional 
retrieval algorithms. Namely, the radiative effect of cloud sides viewed at oblique angles 
seems to contribute most to the enhanced x-values. The results presented here will help 
understand cloud retrieval uncertainties related to cloud inhomogeneity. They 
complement the uncertainty estimates that will start accompanying MODIS cloud 
products in Collection 5 and may eventually help correct for the observed view angle 
dependent biases. 
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Introduction 


Satellite remote sensing is such a complex task that until now it could be done 
only by using one-dimensional (ID) radiative transfer theory, which assumes that cloudy 
pixels are fully covered by horizontally homogeneous clouds and that the pixels’ 
radiative properties are not affected by cloud variability in nearby areas. The use of ID 
theory is often referred to as the plane-parallel approximation. It is true that some 
recently proposed methods (e.g., Marshak et al. 1998, Oreopoulos et al. 2000a, Faure et 
al. 2001, Varnai and Marshak 2002a, Iwabuchi and Hayasaka 2003, Cornet et al. 2004, 
2005) use some aspects of three-dimensional (3D) radiative transfer theory for retrievals 
of cloud optical thickness, but these novel methods are not yet ready for operational use. 

In recent years, several observational studies examined whether ID radiative 
theory gives accurate results in satellite remote sensing. These studies found that, under 
certain conditions, 3D effects cause significant problems. Specifically, they revealed that 
3D effects can make clouds appear too smooth (e.g., Marshak et al. 1995, Davis et al. 
1997), too bright and thick (e.g., Loeb and Davies 1996, Loeb and Coakley 1998), and 
artificially asymmetric (Varnai and Marshak 2002a, b). 

While the papers above focused mainly on overhead satellite views, some studies 
examined 3D effects for oblique views. A comparison of GOES and Meteosat radiances 
for scenes that were viewed from different directions by the two satellites did not reveal 
any influence of 3D effects (Rossow 1989). Using multiangle MISR (Multiangle Imaging 
SpectroRadiometer) observations, however, Horvath and Davies (2004) showed that the 
angular pattern of cloud reflection rarely fits the expectations based on the plane-parallel 
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approximation. Examining ERBE (Earth Radiation Budget Experiment), AVHRR 
(Advanced Very High-Resolution Radiometer), and POLDER (Polarization and 
Directionality of the Earth’s Reflectances) data, some other studies (Loeb and Davies 
1997; Loeb and Coakley 1998; Buriez et al. 2001) found that for low sun, 3D interactions 
such as shadowing make clouds appear too dark from oblique views facing the sun, and 
that this makes ID retrievals underestimate cloud optical thickness. Theoretical studies 
(e.g., Davies 1984; Kobayashi 1993; Loeb et al. 1998; Szczap et al. 2000; Varnai 2000; 
Chambers et al. 2001; Iwabuchi and Hayasaka 2002) have long suggested that 3D effects 
have an opposite influence for oblique views facing away from the sun— but the 
observations cited above have not confirmed unambiguously the existence of this 
enhanced backscatter from sunlit slopes. Zuidema et al. (2003) found that in highly 
inhomogeneous cumulus congestus clouds, oblique backscatter reflectances observed by 
MISR exceeded 3D radiative transfer calculations based on cloud structure retrieved from 
the MISR nadir camera using the plane-parallel approximation. Recently, Marchand and 
Ackerman (2004) found that stratocumulus reflection in backscatter direction was 
stronger in MISR observations than in ID or even 2D simulations based on a variety of 
ground-based and satellite observations. 

Finally, theoretical studies (e.g., Davies 1984, Breon 1992, Kobayashi 1993, 
Iwabuchi and Hayasaka 2002) also showed that cloud inhomogeneities can enhance 
reflection into oblique side scatter directions relative to reflection into nadir direction. 
The observations of Minnis (1989) revealed that cloud side viewing must occur 
frequently, because it increases cloud coverage significantly for oblique views. Still, 
while the viewing of cloud sides can be expected to yield larger retrieved optical 
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thicknesses for oblique side scatter views than for overhead views, the significance of this 
effect has not yet been determined through observations. 

The goal of this paper is to analyze the view angle dependence of a one year long 
MODIS cloud optical thickness dataset, and to examine whether cloud inhomogeneity 
has a significant influence on this view angle dependence. Section 2 describes the data 
we analyzed, Section 3 outlines our methodology, and Section 4 presents the results of 
our analysis. Section 5 then discusses potential reasons for the observed view angle 
dependence, including the effects of cloud sides on the retrievals. Finally, Section 6 
offers a brief summary and discusses the results’ main implications. 

2. Observations 

This study took advantage of the unprecedented abundance of high-quality, easy- 
to-use, and freely available cloud products from new Earth Observing System (EOS) 
satellites. In particular, it used a dataset extracted from the continuous stream of 
incoming MODIS observations at the Goddard Earth Sciences Data and Information 
Services Center (GES DISC) MODIS data pool. The dataset includes observations from 
virtually all daytime granules from the MODIS instruments on both the Terra and Aqua 
satellites for a one year long period ranging from August 2004 to July 2005. The dataset 
includes 1 km-resolution Collection 4 MODIS products such as the 1 1 /<m brightness 
temperature, cloud phase, cloud optical thickness, cloud particle size, and cloud top 
pressure, as well as geolocation parameters such as latitude, longitude, surface type, and 
sun-view geometry. To reduce data volume, these parameters were saved only for about 
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every 14 th row in the MODIS images. To help examine the influence of local cloud 
variability, 1 1 pim brightness temperature and cloud optical thickness values were also 
saved for both neighbors of each row. Finally, we note that in order to avoid the effects of 
uncertainties in cloud detection and in ice crystal scattering phase functions, this study 
analyzed only liquid phase pixels that were flagged as “high confidence” by the 
operational MODIS optical thickness retrieval algorithm. 


3. Methodology 


MODIS is suitable for analyzing the view angle dependence of retrieved cloud 
parameters because clouds are viewed from nadir direction at the MODIS swath center 
and from highly oblique directions at the swath edges, with maximum viewing zenith 
angles exceeding 60°. It is important to note that the oblique views are not aligned with 
the solar azimuth and represent side scattering at both swath edges: At low solar 
elevations, observations are typically from 60° and 110° relative azimuths at the two 
swath edges (Figure 1). 

One approach to identifying the influence of 3D effects is to contrast the view 
angle dependence of optical thicknesses (x) retrieved for homogeneous and 
inhomogeneous clouds, for which 3D effects are expected to be weaker and stronger, 
respectively. Following Varnai and Marshak (2002a), we characterize the degree of 
inhomogeneity at a given pixel through AT, the 1 1 /rm brightness temperature gradient in 
a direction close to the solar azimuth: 


AT - 


T f -T b 


d 


( 1 ) 
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where T is brightness temperature, the subscripts / and b identify the neighboring pixels 
in front and behind our pixel, as viewed from the solar direction, and d is the distance 
between these two neighboring pixels (d = 2 km or, if the solar azimuth is close to 
diagonal in the MODIS image, d=-yj2- 2 km). We note that using the range of brightness 
temperatures in a 3 X 3 pixel window for characterizing cloud variability produced nearly 
identical results to our approach, whereas using x-variability proved less effective in 
separating homogeneous and inhomogeneous clouds for the purpose of creating two 
cloud categories with distinct view angle dependencies, perhaps because optical thickness 
itself is a product of the plane-parallel approximation. 

Using the AT values defined in Eq. (1), we separated cloudy pixels that had high- 
confidence liquid phase x-retrievals into three equally populous categories based on the 
degree of local cloud variability. Cloudy pixels over ocean and land were assigned into 
the most homogeneous category if their AT value was less than 0.2 °C/km or 0.5 °C/km, 
respectively. Cloudy pixels over ocean and land were assigned into the most 
inhomogeneous category if their AT value was greater than 0.9 °C/km or 2.2 °C/km, 
respectively. In summary: 

For ocean: 

Homogeneous: AT < 0.2 °C/km 
Intermediate: 0.2 °C/km < AT < 0.9 °C/km 
Inhomogeneous: AT > 0.9 °C/km 
For land: 

Homogeneous: AT < 0.5 °C/km 
Intermediate: 0.5 °C/km < AT < 2.2 °C/km 
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Inhomogeneous: A T> 2.2 °C/km 

Local temperature gradients are higher over land than over ocean probably 
because stronger surface heating causes stronger convection over land. Finally, we note 
that in plots that do not separate clouds over land and ocean, we combine the two 
categories by weighting them according to the number of pixels they contain. 

4. Results 

Considering a variety of solar zenith angles (9 0 ), Figure 2 shows the way the mean 
retrieved x-values change with view angle (0) for the most homogeneous and most 
inhomogeneous third of cloudy pixels. The figure clearly confirms the findings of Loeb 
and Davies (1996) and Loeb and Coakley (1998) in that 3D effects cause the retrieved x- 
values to increase with solar zenith angle. It also shows that homogeneous clouds tend to 
be thicker than inhomogeneous clouds. In addition. Panel a indicates that for 
homogeneous clouds, the plane-parallel approximation produces consistent results that 
don’t change much with view direction— although both homogeneous and 
inhomogeneous clouds tend to be slightly thicker at the more forward scattering swath 
edge than at the opposite edge. While we have not examined the reasons for this slight 
cross-track trend, it appears possible that the trend reflects true changes in cloud 
properties, which could arise from a combination of systematic cross-track variations in 
local time and latitude (due to the sun-synchronous orbits of the Terra and Aqua 
satellites) on one hand, and the latitudinal distribution of cloud properties and the daily 
cycle of cloud development on the other. Panel b, however, reveals that for 
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inhomogeneous clouds, the plane-parallel approximation yields substantially higher x- 
values for oblique views than for overhead views if the sun is fairly oblique. Note that for 
the most oblique solar zenith angles, the difference exceeds 50%. 

Figure 3 shows that clouds tend to be thicker over land than over ocean, and also 
that the U-shape for inhomogeneous clouds is more pronounced over land than over 
ocean. This is consistent with the AT local temperature gradients being larger over land 
than over ocean, perhaps because stronger surface heating causes stronger convection 
over land. This tendency, however, differs from the results of Oreopoulos and Cahalan 
(2005), who found that large-scale variability over 1° by 1° areas tends to be stronger 
over ocean. The opposite tendencies at small and large scales indicate that while cloud 
fields tend to be bumpier over land than over ocean, their statistical properties vary more 
gradually over land. 

Figure 4 illustrates the way the depth of the U-shape changes with solar zenith 
angle for clouds over land and ocean. The figure characterizes the depth of the U-shape 
through the c(0 o ) coefficient obtained by fitting to the mean optical thicknesses in 5°- 
wide 0 O intervals a second order polynomial in the form of: 

f(e,d 0 ) = a{d 0 ) + b{d 0 ) • 6+ c(e 0 y e 2 . (2) 

Figure 5 indicates that the U-shape is stronger for high water clouds, which tend 
to have larger variability than low clouds. This is consistent with high clouds being 
thicker both optically and geometrically, which allows more pronounced 
inhomogeneities. Figure 6 shows that the difference between optical thicknesses retrieved 
at overhead and oblique views is significant throughout the entire range of cloud 
thicknesses: for oblique observations, optical thicknesses smaller and larger than 1 1 are 
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less and more frequent, respectively. Finally, the results (not shown) indicate that the U- 
shape is similar over the Northern and Southern hemispheres, for the Terra and Aqua 
satellites, and for various seasons throughout the year. This allows us to conclude that the 
presence of the U-shape is not restricted to a particular cloud type or location, but 
represents the general behavior of MODIS retrievals for inhomogeneous clouds. 

5. Potential reasons for the observed behaviors 

This section examines the following potential explanations for the U-shape in 
Figure 3: 

• whether inhomogeneous clouds viewed obliquely are indeed thicker than those 
viewed from overhead, 

• whether inhomogeneous clouds behave differently than homogeneous clouds 
because of their different altitude, 

• whether uncertainties in surface reflection, cloud phase, or cloud altitude may 
cause the behaviors in Figure 3, 

• whether cross-track changes in MODIS pixel size may explain the observations, 

• whether horizontal photon transport in bumpy clouds may cause the U-shapes in 
Figure 3, 

• and finally, whether the viewing of cloud sides may cause the U-shapes. 

We first examine whether the U-shape in Figure 3 may reflect the true behavior of 
inhomogeneous clouds. It appears unlikely that cross-track variations in local time or 
latitude— combined with the latitudinal distribution of cloud properties and the daily 
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cycle of cloud development— could cause the U-shapes observed for inhomogeneous 
clouds, because the curves remain very similar for various combinations of satellite, 
hemisphere, and season, even though the local times and latitudes of observations are 
quite different for the various combinations. Thus it is more likely that the U-shapes do 
not reflect the true behavior of inhomogeneous clouds, but are caused by some artifact in 
the retrievals instead. 

The next question to consider is whether the difference between the behaviors of 
homogeneous and inhomogeneous clouds is caused by inhomogeneity itself or by some 
other difference between the populations of homogeneous and inhomogeneous clouds. 
Figures 3 and 5 imply that the U-shapes do not arise from inhomogeneous clouds 
occurring over different surfaces or at different altitudes than homogeneous clouds. 
Another possibility is that if, due to different updraft speeds, ice crystals had different 
shapes in homogeneous and inhomogeneous clouds, or if cloud inhomogeneity made it 
more difficult to detect cloud phase accurately, ice contamination in our supposedly 
liquid cloud dataset could cause different view-angle dependencies for homogeneous and 
inhomogeneous clouds. However, this is also unlikely since, as shown in Figure 7, the U- 
shape is present even if only warm pixels with brightness temperatures exceeding 0°C 
are considered. 

As a result of the considerations above, it appears very likely that the U-shapes in 
Figure 3 are indeed caused by cloud inhomogeneity. One potential mechanism for this 
would be if cloud top pressure retrievals were influenced by cloud variability over the 
operational retrievals’ 5 km by 5 km domain, and this caused errors in atmospheric 
correction over inhomogeneous clouds. Because absorption by tropospheric gases is 


9 



negligible at visible MODIS wavelengths, the effects would be strongest in the correction 
for Rayleigh scattering. However, simple ID calculations indicate that errors in cloud 
altitude could cause much weaker effects than those observed, and the changes would 
have opposite sign at the two sides of MODIS swaths — which is inconsistent with the 
observations always indicating higher-than-nadir values for oblique views. 

Increases in MODIS pixel size for oblique views can also influence the view- 
angle dependence of retrieved optical thicknesses, because averaging of radiances over 
larger areas cause stronger plane-parallel biases (e.g., Oreopoulos and Davies 1998). 
Because of the concavity of the ID reflectance vs. x curve, however, averaging always 
decreases the retrieved x- values, and so stronger averaging at oblique views would create 
a D shape rather than the U-shape observed in Figure 3. 

Theoretical simulations (e.g., Davies 1984, Breon 1992, Kobayashi 1993) indicate 
that horizontal photon transport can result in larger optical thicknesses being retrieved for 
the oblique views typical of MODIS observation geometry than for overhead views, but it 
is unclear whether this tendency would disappear for high sun. 3D effects could also 
cause the enhancement at oblique views in Figure 3 by increasing the variability of the 
reflectance field: Because of the nonlinearity of the ID reflectance vs. x curve, retrieved 
x-values change more if 3D effects increase, rather than decrease reflectance by a certain 
amount. As a result, retrieved x-values can increase more on sunlit slopes than they 
decrease in shadowy slopes, and this can increase the average x- value even if 3D effects 
did not increase the average reflectance. While this behavior has been reproduced in our 
3D Monte Carlo simulations, the increase for oblique views occurred exclusively at 
pixels whose brightness was enhanced by 3D effects. This, however, does not seem to fit 
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the observations in which the U-shape appears even for pixels for which T { < T b , which 
tend to lie on slopes facing away from the sun and to be darkened by 3D effects (Fig. 8). 

Finally, cloud side viewing also influences the view angle dependence of 
retrieved x- values. The observations of Minnis (1989) indicate that cloud side viewing 
must occur quite frequently indeed, because it significantly increases cloud coverage for 
oblique views. Reflection into oblique directions is generally enhanced through cloud 
sides, including sides of thicker cloud elements that are surrounded by significantly 
thinner regions. Cloud side viewing can also cause the U-shape in Figure 3 because while 
overhead views can see through small subpixel gaps in cloudiness, oblique'views tend to 
see cloud sides instead. This can result in stronger plane-parallel biases and thus stronger 
reductions in retrieved optical thicknesses for overhead views than for oblique views. 
While such side-viewing occurs for any solar elevation, its effects are stronger for 
oblique sun because of the more pronounced the nonlinearity (i.e., earlier saturation) of 
the ID reflectance vs. x curve. Still, we are somewhat cautioned by the U-shape being 
present even for pixels that are surrounded by cloudy pixels on all sides and are farther 
than 10 km away from clear pixels in cross-track direction (Figure 9). This observation, 
however, doesn’t exclude side viewing from being an important factor in creating the U- 
shape in Figure 3, because subpixel gaps can occur in areas that appear overcast at 1 km 
resolution, and also because side viewing can have important effects in multilayer cloud 
systems or even in highly variable clouds. 
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6. Summary and discussion 


This paper examined whether cloud inhomogeneity introduces any view angle 
dependent biases into MODIS cloud optical thickness (x) retrievals, which use the plane- 
parallel approximation and hence assume cloud homogeneity. The influence of cloud 
inhomogeneity was identified by contrasting the view angle dependence of mean x values 
retrieved for clouds that were deemed homogeneous or inhomogeneous based on the 
local gradient in 1 1 pim brightness temperature. 

The analysis of liquid phase clouds in a one year-long global dataset of 
Collection 4 MODIS cloud products revealed that while optical thickness retrievals give 
remarkably consistent results at all view directions for homogeneous clouds, they give 
systematically higher x-values at oblique views than at overhead views for 
inhomogeneous clouds if the sun is fairly oblique. The mean optical thickness retrieved 
for the most inhomogeneous third of cloudy pixels is more than 30% higher for oblique 
views at swath edges than for overhead views at the swath center if the solar zenith angle 
is 60°; the difference exceeds 40% if the solar zenith angle is 70°. The observations 
reveal that the dependence on view angle is stronger for higher clouds and for clouds 
over land, that it is present over a wide range of cloud thicknesses at both hemispheres 
through all seasons, and that it is quite similar in observations by the Terra and Aqua 
satellites. 

After considering several potential scenarios, the paper concluded that the 
observed behavior is indeed caused by cloud inhomogeneities influencing ID cloud 
property retrievals, and not by other differences between homogeneous and 
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inhomogeneous clouds that are unrelated to inhomogeneity itself (e.g., in microphysics). 
The paper discussed several mechanisms through which cloud inhomogeneity may 
influence the view angle dependence of x-values. The most likely candidate is the 
increased viewing of cloud sides from oblique directions. Cloud sides can enhance 
reflection into oblique directions, and while cloudy pixels may contain small gaps among 
broken clouds in overhead views, these dark gaps tend to be filled by cloud sides in 
oblique views. 

Once the dominant mechanism is identified unambiguously, it will be possible to 
determine whether the view angle dependent biases decrease x-values for nadir views or 
increase them for oblique views. We note, however, that the view angle dependent biases 
discussed in this paper are only one component of the overall radiative effect of cloud 
inhomogeneity which, as discussed in earlier studies (Loeb and Davies 1996; Loeb and 
Coakley 1998), increases x with solar zenith angle even for overhead views— though by 
not as much as for oblique views. 

The results presented here can help improve future versions of the x-retrieval 
uncertainty estimates that will start accompanying MODIS cloud products in Collection 
5. These uncertainty estimates consider only factors within a ID framework (such as 
uncertainties in calibration, in atmospheric correction, and in surface albedo), whereas 
our results suggest that identifying inhomogeneous pixels through local brightness 
temperature gradients could help incorporating view angle dependent cloud 
inhomogeneity effects as well— and this approach may eventually help correct for the 
observed view angle dependent biases. 
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Figure 1. Schematic view of MODIS observational geometry. 
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Figure 2. View angle-dependence of mean retrieved optical thickness. Only liquid phase 
clouds with high-confidence retrievals are considered. Each line represents a separate 
solar zenith angle (0 O ) interval, (a) Most homogeneous third of cloudy pixels; (b) Most 
inhomogeneous third of cloudy pixels. 
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Figure 3. View angle dependence of mean optical thickness of inhomogeneous clouds 
over land and ocean. For increased clarity, each curve represents the average for 5 
different solar zenith angle intervals ranging from 55° to 80°. The thick solid lines 
represent second-order polynomial fits to the curves. 
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obtained by fitting second-order polynomial to the view-angle dependence of retrieved x- 
values. 
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Figure 5. View angle dependence of mean x for the most homogeneous and most 
inhomogeneous third of cloudy pixels that have MODIS-estimated cloud top pressures 
below and above 700 hPa. Accordingly, the “high” and “low” curves represent clouds 
with tops higher and lower than about 3 km, respectively. For increased clarity, each 
curve represents the average for 5 different solar zenith angle intervals ranging from 55° 
to 80°. The figure also displays the c coefficients obtained by fitting a second-order 
polynomial to each curve. 
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Figure 6. Probability distribution function (PDF) of x for the most inhomogeneous third 


of cloudy pixels for nadir view and for oblique views slightly oriented toward forward 


and back scatter. The viewing zenith angle is in the 50° -60° and the 0°-5° range, 


respectively; the solar zenith angle is between 60° and 70°. 


25 





Back scatter 


Viewing zenith angle (°) 

Forward scatter 


Figure 7. View angle dependence of mean x for inhomogeneous pixels with 1 1 pim 
brightness temperatures warmer than 0°C, for solar zenith angle ranging from 55° to 80°. 
The thick line represents a second order polynomial fit to the data, with c- 1.0* 10' 3 . 
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Figure 8. View angle dependence of x for inhomogeneous pixels that lie on slopes tilted 
toward and away from the sun ( T f > T b and T { < T b , respectively). The curves represent the 
average for 5 different solar zenith angle intervals ranging from 55° to 80°. The quadratic 
polynomial coefficients for slopes tilted toward and away from the sun are c=1.82*10' 3 
and 1.28* 10' 3 , respectively. 
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Figure 9. View angle dependence of mean optical thickness for inhomogeneous pixels 
that occur at various distances D from the nearest cloud-free pixel in cross- track 
direction. The curves represent the average for 5 different solar zenith angle intervals 
ranging from 55° to 80°. Values of quadratic polynomial coefficient c are also displayed. 
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